Sex differences orchestrated by androgens at single-cell resolution.

Sex differences in mammalian complex traits are prevalent and are intimately associated with androgens1-7. However, a molecular and cellular profile of sex differences and their modulation by androgens is still lacking. Here we constructed a high-dimensional single-cell transcriptomic atlas comprising over 2.3 million cells from 17 tissues in Mus musculus and explored the effects of sex and androgens on the molecular programs and cellular populations. In particular, we found that sex-biased immune gene expression and immune cell populations, such as group 2 innate lymphoid cells, were modulated by androgens. Integration with the UK Biobank dataset revealed potential cellular targets and risk gene enrichment in antigen presentation for sex-biased diseases. This study lays the groundwork for understanding the sex differences orchestrated by androgens and provides important evidence for targeting the androgen pathway as a broad therapeutic strategy for sex-biased diseases.

Mutant U2AF1-induced mis-splicing of mRNA translation genes confers resistance to chemotherapy in acute myeloid leukemia.

Patients with primary refractory acute myeloid leukemia (AML) have a dismal long-term prognosis. Elucidating the resistance mechanisms to induction chemotherapy could help identify strategies to improve AML patient outcomes. Herein, we retrospectively analyzed the multi-omics data of more than 1,500 AML cases and found that patients with spliceosome mutations had a higher risk of developing refractory disease. RNA splicing analysis revealed that the mis-spliced genes in refractory patients converged on translation-associated pathways, promoted mainly by U2AF1 mutations. Integrative analyses of binding and splicing in AML cell lines substantiated that the splicing perturbations of mRNA translation genes originated from both the loss and gain of mutant U2AF1 binding. In particular, the U2AF1-S34F and U2AF1-Q157R mutants orchestrated the inclusion of exon 11 (encoding a premature termination codon) in the eukaryotic translation initiation factor 4A2 (EIF4A2). This aberrant inclusion led to reduced eIF4A2 protein expression via nonsense-mediated mRNA decay. Consequently, U2AF1 mutations caused a net decrease in global mRNA translation that induced the integrated stress response (ISR) in AML cells, which was confirmed by single-cell RNA-seq. The induction of ISR enhanced the ability of AML cells to respond and adapt to stress, contributing to chemoresistance. A pharmacologic inhibitor of ISR, ISRIB, sensitized U2AF1 mutant cells to chemotherapy. These findings highlight a resistance mechanism by which U2AF1 mutations drive chemoresistance and provide a therapeutic approach for AML through targeting the ISR pathway.

The prognostic value of CD8+ CTLs, CD163+ TAMs, and PDL1 expression in the tumor microenvironment of primary central nervous system lymphoma.

To explore immune cell infiltration and PDL1 expression in the tumor microenvironment (TME) of primary central nervous system lymphoma (PCNSL), we performed immunohistochemical staining on paraffin-embedded tumor tissues from 34 patients diagnosed with PCNSL. CD8 and CD163 positive cells were manually counted, and PDL1 expression was quantified by the H-score scoring method in the tumor center and around the tumor. The Kaplan-Meier method was used to analyze the prognostic value of the TME. We found obvious infiltration of CD8+ CTLs and CD163+ TAMs in the TME of PCNSL patients. And PDL1 was expressed in the tumor center as well as around the tumor. Survival analysis showed that high CD8+ CTLs levels and high intratumoral PDL1 expression were significantly correlated with longer OS. High CD8+ CTLs and CD163+ TAMs levels were associated with longer PFS.

A prediction model of elderly hip fracture mortality including preoperative red cell distribution width constructed based on the random survival forest (RSF) and Cox risk ratio regression.

An independent correlation between pre-RDW and 1-year mortality after surgery in elderly hip fracture can be used to predict mortality in elderly hip fracture patients and has predictive significance in anemia patients. With further research, a treatment algorithm can be developed to potentially identify patients at high risk of preoperative mortality. INTRODUCTION: Red blood cell distribution width (RDW) is an independent predictor of various disease states in elderly individuals, but its association with the prognosis of elderly hip fracture patients is controversial. This study aimed to evaluate the prognostic value of RDW in such patients, construct a prediction model containing RDW using random survival forest (RSF) and Cox regression analysis, and compare RDW in patients with and without anemia. METHODS: We retrospectively analyzed the data of elderly patients who underwent hip fracture surgery, selected the best variables using RSF, stratified the independent variables by Cox regression analysis, constructed a 1-year mortality prediction model of elderly hip fracture with RDW, and conducted internal validation and external validation. RESULTS: Two thousand one hundred six patients were included in this study. The RSF algorithm selects 12 important influencing factors, and Cox regression analysis showed that eight variables including preoperative RDW (pre-RDW) were independent risk factors for death within 1-year after hip fracture surgery in elderly patients. Stratified analysis showed that pre-RDW was still independently associated with 1-year mortality in the non-anemia group and not in the anemia group. The nomogram prediction model had high differentiation and fit, and the prediction model constructed by the total cohort of patients was also used for validation of patients in the anemia patients and obtained good clinical benefits. CONCLUSION: An independent correlation between pre-RDW and 1-year mortality after surgery in elderly hip fracture can be used to predict mortality in elderly hip fracture patients and has predictive significance in anemia patients.

Role and application of chemokine CXCL13 in central nervous system lymphoma.

Chemokine ligand 13 (CXCL13) and its chemokine receptor 5 (CXCR5) both play significant roles in the tumor microenvironment (TME). CXCL13 in cerebrospinal fluid (CSF) has recently been found to have significant diagnostic and prognostic value in primary and secondary central nervous system (CNS) diffuse large B-cell lymphoma (DLBCL), and the CXCL13-CXCR5 axis has been shown to play an important chemotactic role in the TME of CNS-DLBCL. In this review, we first describe the clinical value of CXCL13 in CSF as a prognostic and diagnostic biomarker for CNS-DLBCL. In addition, this review also discusses the specific mechanisms associated with the CXCL13-CXCR5 axis in tumor immunity of primary diffuse large B cell lymphoma of the central nervous system (PCNS-DLBCL) by reviewing the specific mechanisms of this axis in the immune microenvironment of DLBCL and CNS inflammation, as well as the prospects for the use of CXCL13-CXCR5 axis in immunotherapy in PCNS-DLBCL.

[The Prognostic Value of Prognostic Nutritional Index Combined with D-dimer in Patients with Diffuse Large B-Cell Lymphoma].

OBJECTIVE: To explore the effects of prognostic nutritional index (PNI) combined with D-dimer on the prognosis of patients with newly diagnosed diffuse large B-cell lymphoma (DLBCL). METHODS: The clinical data of 73 DLBCL patients at initial diagnosis were retrospectively evaluated, and the optimal cut-off point of PNI and D-dimer were determined by ROC curve. The overall survival (OS) rate and progression-free survival (PFS) rate in different subgroups were compared using Kaplan-Meier survival curves. Univariate and multivariate Cox regression analysis was performed to identify the factors associated with OS. RESULTS: Compared with the low PNI group (PNI<44.775), the high PNI group (PNI≥44.775) had better OS (P =0.022) and PFS (P =0.029), the 2-year OS rates of the two groups were 55.6% and 78.3% respectively (P =0.041). Compared with the high D-dimer group (D-dimer≥0.835), the low D-dimer group (D-dimer<0.835) had better OS (P <0.001) and PFS (P <0.001), the 2-year OS rates of the two groups were 51.4% and 86.8% respectively (P =0.001). Meanwhile, patients in the high PNI+ low D-dimer group had better OS (P =0.003) and PFS (P <0.001) than the other three groups, the 2-year OS rate was statistically different from the other three groups (P <0.05). The multivariate analysis revealed that NCCN-IPI (HR =2.083, 95%CI : 1.034-4.196, P =0.040), PNI (HR =0.267, 95%CI : 0.076-0.940, P =0.040) and PNI+D-dimer (HR =9.082, 95%CI : 1.329-62.079, P =0.024) were the independent risk factors affecting OS in patients with DLBCL. Subgroup analysis showed that PNI, D-dimer, and PNI combined with D-dimer could improve the prognostic stratification in low and low-intermediate risk DLBCL patients. CONCLUSION: High PNI, low D-dimer and combination of high PNI and low D-dimer at initial diagnosis suggest a better prognosis in DLBCL patients.

Identification and Characterization of CD8+ CD27+ CXCR3- T Cell Dysregulation and Progression-Associated Biomarkers in Systemic Lupus Erythematosus.

Systemic Lupus Erythematosus (SLE) etiopathogenesis highlights the contributions of overproduction of CD4+ T cells and loss of immune tolerance. However, the involvement of CD8+ T cells in SLE pathology and disease progression remains unclear. Here, the comprehensive immune cell dysregulation in total 263 clinical peripheral blood samples composed of active SLE (aSLE), remission SLE (rSLE) and healthy controls (HCs) is investigated via mass cytometry, flow cytometry and single-cell RNA sequencing. This is observed that CD8+ CD27+ CXCR3- T cells are increased in rSLE compare to aSLE. Meanwhile, the effector function of CD8+ CD27+ CXCR3- T cells are overactive in aSLE compare to HCs and rSLE, and are positively associated with clinical SLE activity. In addition, the response of peripheral blood mononuclear cells (PBMCs) is monitored to interleukin-2 stimulation in aSLE and rSLE to construct dynamic network biomarker (DNB) model. It is demonstrated that DNB score-related parameters can faithfully predict the remission of aSLE and the flares of rSLE. The abundance and functional dysregulation of CD8+ CD27+ CXCR3- T cells can be potential biomarkers for SLE prognosis and concomitant diagnosis. The DNB score with accurate prediction to SLE disease progression can provide clinical treatment suggestions especially for drug dosage determination.

CD39 (ENTPD1) in tumors: a potential therapeutic target and prognostic biomarker.

As a regulator of the dynamic balance between immune-activated extracellular ATP and immunosuppressive adenosine, CD39 ectonucleotidase impairs the ability of immune cells to exert anticancer immunity and plays an important role in the immune escape of tumor cells within the tumor microenvironment. In addition, CD39 has been studied in cancer patients to evaluate the prognosis, the efficacy of immunotherapy (e.g., PD-1 blockade) and the prediction of recurrence. This article reviews the importance of CD39 in tumor immunology, summarizes the preclinical evidence on targeting CD39 to treat tumors and focuses on the potential of CD39 as a biomarker to evaluate the prognosis and the response to immune checkpoint inhibitors in tumors.

HLA-I Evolutionary Divergence Confers Response to PD-1 Blockade plus Chemotherapy in Untreated Advanced Non-Small Cell Lung Cancer.

PURPOSE: PD-1 blockade plus chemotherapy has become the new standard of care in patients with untreated advanced non-small cell lung cancer (NSCLC), whereas predictive biomarkers remain undetermined. EXPERIMENTAL DESIGN: We integrated clinical, genomic, and survival data of 427 NSCLC patients treated with first-line PD-1 blockade plus chemotherapy or chemotherapy from two phase III trials (CameL and CameL-sq) and investigated the predictive and prognostic value of HLA class I evolutionary divergence (HED). RESULTS: High HED could predict significantly improved objective response rate (ORR), progression-free survival (PFS), and overall survival (OS) in those who received PD-1 blockade plus chemotherapy [in the CameL trial, ORR: 81.8% vs. 53.2%; P = 0.032; PFS: hazard ratio (HR), 0.47; P = 0.012; OS: HR, 0.40; P = 0.014; in the CameL-sq trial, ORR: 89.2% vs. 62.3%; P = 0.007; PFS: HR, 0.49; P = 0.005; OS: HR, 0.38; P = 0.002], but not chemotherapy. In multivariate analysis adjusted for PD-L1 expression and tumor mutation burden, high HED was independently associated with markedly better ORR, PFS, and OS in both trials. Moreover, the joint utility of HED and PD-L1 expression showed better performance than either alone in predicting treatment benefit from PD-1 blockade plus chemotherapy. Single-cell RNA sequencing of 58,977 cells collected from 11 patients revealed that tumors with high HED had improved antigen presentation and T cell-mediated antitumor immunity, indicating an inflamed tumor microenvironment phenotype. CONCLUSIONS: These findings suggest that high HED could portend survival benefit in advanced NSCLC treated with first-line PD-1 blockade plus chemotherapy. See related commentary by Dimou, p. 4706.

Comparative Analysis of the Quality in Ripe Fruits of Cuiguan Pear from Different Regions.

The Cuiguan pear is called "June snow" and the skin is thin; the meat is crisp and juicy; the taste is thick and fresh; and the juice is rich and sweet. In this study, the volatile organic compounds and the sensory and physicochemical parameters of the Cuiguan pear from four different regions of China (Sichuan (SC), Shangdong (SD), Chongming (CM), Zhuanghang (ZH)) were assessed. The highest differences in the physicochemical parameters were observed between four regions. The volatile fingerprints of GC-IMS showed great differences in the volatile of the Cuiguan pear, which suggested that the aroma of pears could be largely impacted by origin areas. (E)-ethyl-2-hexenoate can be used to distinguish between the 'CM' and pears from other regions. High contents of 2-heptanone, 1-pentanol, 1-butanol, 3-methylbutanol, butyl 2-methylbutanoate, heptyl acetate and butyl acetate were observed in the 'SD'. Dimethyl trisulfide, 6-methyl-5-hepten-2-one, 3-hydroxy-2-butanone, 1-penten-3-one, beta-pinene, γ-terpinene, propanal, (e)-2-pentenal, (e)-2-heptenal, 1-pentanol and 3-methyl-1-pentanol were primarily contained in the 'ZH'. Principal component analysis showed that there was very good discrimination based on the information obtained from GC-IMS for four samples. These findings were in agreement with the sensory analysis. In the opinion of the respondents to the consumer test, 'ZH' resulted in the most appreciated sample based on the average scores of the acceptability. This study provides some reference for the development and utilization of the Cuiguan pear.

Frequent Gene Mutations and Their Possible Roles in the Pathogenesis, Treatment and Prognosis of Primary Central Nervous System Lymphoma.

Primary central nervous system lymphoma (PCNSL) is a rare extranodal non-Hodgkin lymphoma with poor prognosis. In recent years, the emergence of genetic subtypes of systematic diffuse large B-cell lymphoma has highlighted the importance of molecular genetics, but large-scale research on the molecular genetics of PCNSL is lacking. Herein, we summarize the frequent gene mutations and discuss the possible pathogenesis of PCNSL. Myeloid differentiation primary response gene 88 (MYD88) and CD79B mutations, which cause abnormal activation of noncanonical nuclear factor-κB, are prominent genetic abnormalities in PCNSL. They are considered to play a major role in the pathogenesis of PCNSL. Other genes, such as caspase recruitment domain family member 11 (CARD11), tumor necrosis factor alpha induced protein 3 (TNFAIP3), transducin (ß)-like 1 X-linked receptor 1, cyclin dependent kinase inhibitor 2A, PR domain zinc finger protein 1, and proviral insertion in murine malignancies 1, are also frequently mutated. Notably, the pathogenesis of immune insufficiency-associated PCNSL is related to Epstein-Barr virus infection, and its progression may be affected by different signaling pathways. The different mutational patterns in different studies highlight the heterogeneity of PCNSL. However, existing research on the molecular genetics of PCNSL is still limited, and further research into PCNSL is required to clarify the genetic characteristics of PCNSL.

[Prognostic Value of A New Risk Stratification Model (R2-ISS) in Newly Diagnosed Multiple Myeloma].

OBJECTIVE: To investigate the significance of a new risk stratification model (R2-ISS) in evaluating the prognosis of newly diagnosed multiple myeloma (MM). METHODS: Clinical data of 116 newly diagnosed MM patients admitted to Lanzhou University Second Hospital from June 2012 to March 2021 were retrospectively analyzed. According to R2-ISS, these patients were divided into four groups: low risk, low-intermediate risk, intermediate-high risk, and high risk. The significance of R2-ISS on prognosis of the patients was analyzed. RESULTS: Survival analysis showed that R2-ISS was associated with progression-free survival (PFS) (P=0.042) and overall survival (OS) (P=0.014). Cox univariate analysis showed that lactate dehydrogenase, serum calcium, serum creatinine, ß2-microglobulin, ISS, R-ISS, R2-ISS, t(4;14), and autologous hematopoietic stem cell transplantation (ASCT) were the influencing factors of OS in newly diagnosed MM patients (all P<0.05). Cox multivariate analysis showed that R-ISS, R2-ISS, and ASCT were independent risk factors affecting OS (all P<0.05). In addition, survival analysis of patients with different R2-ISS showed that ASCT improved PFS and OS. CONCLUSION: R2-ISS has prognostic value for newly diagnosed MM patients, while ASCT can improve the prognosis of patients with different R2-ISS.

Feasibility of Left Bundle Branch Area Pacing Combined with Atrioventricular Node Ablation in Atrial Fibrillation Patients with Heart Failure.

BACKGROUND: Pacemaker implantation combined with atrioventricular node ablation (AVNA) could be a practical choice for atrial fibrillation (AF) patients with heart failure (HF). Left bundle branch area pacing (LBBaP) has been widely reported. OBJECTIVES: To explore the safety and efficacy of LBBaP combined with AVNA in AF patients with HF. METHODS AND RESULTS: Fifty-six AF patients with HF attempted LBBaP and AVNA from January 2019 to December 2020. Standard LBBaP was achieved in forty-six patients, and another ten received left ventricular septal pacing (LVSP). The cardiac function indexes and pacemaker parameters were evaluated at baseline, and we conducted a 1-month and 1-year follow-up. RESULT: At the time of implantation and 1-month and 1-year follow-up, QRS duration of LVSP group was longer than that of LBBaP group. The pacemaker parameters remained stable in both the LBBaP and LVSP groups. At 1-month and 1-year follow-up after LBBaP and AVNA, left ventricular ejection fraction, left ventricular end-diastolic diameter, and NYHA classification continued to improve. Baseline left ventricular ejection fraction and QRS duration change at implantation can predict the magnitude of improvement of left ventricular ejection fraction at 1-year after LBBaP. Baseline right atrial left-right diameter, the degree of tricuspid regurgitation, and interventricular septum thickness may be the factors affecting the success of LBBaP. CONCLUSION: LBBaP combined with AVNA is safe and effective for patients with AF and HF. Baseline right atrial left-right diameter, the degree of tricuspid regurgitation, and interventricular septum thickness may be the factors affecting the success of LBBaP.

Large-Scale In Vitro and In Vivo CRISPR-Cas9 Knockout Screens Identify a 16-Gene Fitness Score for Improved Risk Assessment in Acute Myeloid Leukemia.

PURPOSE: The molecular complexity of acute myeloid leukemia (AML) presents a considerable challenge to implementation of clinical genetic testing for accurate risk stratification. Identification of better biomarkers therefore remains a high priority to enable improving established stratification and guiding risk-adapted therapy decisions. EXPERIMENTAL DESIGN: We systematically integrated and analyzed the genome-wide CRISPR-Cas9 data from more than 1,000 in vitro and in vivo knockout screens to identify the AML-specific fitness genes. A prognostic fitness score was developed using the sparse regression analysis in a training cohort of 618 cases and validated in five publicly available independent cohorts (n = 1,570) and our RJAML cohort (n = 157) with matched RNA sequencing and targeted gene sequencing performed. RESULTS: A total of 280 genes were identified as AML fitness genes and a 16-gene AML fitness (AFG16) score was further generated and displayed highly prognostic power in more than 2,300 patients with AML. The AFG16 score was able to distill downstream consequences of several genetic abnormalities and can substantially improve the European LeukemiaNet classification. The multi-omics data from the RJAML cohort further demonstrated its clinical applicability. Patients with high AFG16 scores had significantly poor response to induction chemotherapy. Ex vivo drug screening indicated that patients with high AFG16 scores were more sensitive to the cell-cycle inhibitors flavopiridol and SNS-032, and exhibited strongly activated cell-cycle signaling. CONCLUSIONS: Our findings demonstrated the utility of the AFG16 score as a powerful tool for better risk stratification and selecting patients most likely to benefit from chemotherapy and alternative experimental therapies.

Single-cell entropy network detects the activity of immune cells based on ribosomal protein genes.

We developed a new computational method, Single-Cell Entropy Network (SCEN) to analyze single-cell RNA-seq data, which used the information of gene-gene associations to discover new heterogeneity of immune cells as well as identify existing cell types. Based on SCEN, we defined association-entropy (AE) for each cell and each gene through single-cell gene co-expression networks to measure the strength of association between each gene and all other genes at a single-cell resolution. Analyses of public datasets indicated that the AE of ribosomal protein genes (RP genes) varied greatly even in the same cell type of immune cells and the average AE of RP genes of immune cells in each person was significantly associated with the healthy/disease state of this person. Based on existing research and theory, we inferred that the AE of RP genes represented the heterogeneity of ribosomes and reflected the activity of immune cells. We believe SCEN can provide more biological insights into the heterogeneity and diversity of immune cells, especially the change of immune cells in the diseases.

One-Year Recording of Cardiac Arrhythmias in a Non-Infected Population with Cardiac Implantable Devices During the COVID-19 Pandemic.

OBJECTIVE: Coronavirus disease 2019 (COVID-19) was associated with a higher risk of arrhythmia in infected patients. However, there are no reports about the effect of the ongoing pandemic on arrhythmias in the non-infected population. We measured the arrhythmia burden in a non-infected population with cardiac implantable devices. METHODS: The arrhythmia burden during the COVID-19 pandemic was compared to a 6-month interval in the pre-COVID-19 period. The COVID-19 pandemic was divided into high-risk (17 January 2020 to 16 March 2020) and low-risk periods (17 March 2020 to 17 July 2020) according to whether there were locally infected patients. Arrhythmia burdens were compared among the pre-COVID-19, high-risk, and low-risk periods. RESULTS: A total of 219 patients with 1859 episodes were included. We observed a larger proportion of patients with atrial fibrillation (AF) during the COVID-19 pandemic (38.36% vs 26.03%, p = 0.006). There was not significantly more ventricular arrhythmia during the COVID period than the pre-COVID-19 period (p > 0.05). During the high-risk period, daily frequency of non-sustained ventricular tachycardia (NSVT) (0.0172, 0.0475 vs 0.0109, 0.0164, p < 0.05), atrial tachycardia (AT) (0.0345, 0.0518 vs 0.0164, 0.0219 p < 0.05) and AF (0.0345, 0.0432 vs 0.0164, 0.0186, p < 0.05) and daily duration of NSVT (0.1982, 0.2845 vs 0.0538, 0.1640 p < 0.05) were higher and longer than those in the pre-COVID-19 period. Regression modeling showed that the impact of COVID-19 pandemic lead to an increased onset of AF (odds ratio 2.465; p < 0.01). Patients with paroxysmal AF who had undergone a previous radiofrequency ablation had a lower burden of AF (incidence 21.43% vs 55.00%, P = 0.049, daily frequency 0.0000, 0.0027 vs 0.0000, 241.7978, P = 0.020) during the pandemic. CONCLUSION: The COVID-19 pandemic contributed to a higher burden of arrhythmias in non-infected patients. Patients would experience a lower burden of AF following radiofrequency ablation treatment, and this effect persisted during the pandemic.

Catheter Ablation of Ventricular Arrhythmias Originating From the Region of DGCV-AIV via a Swartz Sheath Support Approach.

Aims: This study aimed to investigate an appropriate catheter manipulation approach for ventricular arrhythmias (VAs) originating from the left ventricular epicardium adjacent to the transitional area from the great cardiac vein to the anterior interventricular vein (DGCV-AIV). Methods: A total of 123 patients with DGCV-AIV VAs were retrospectively analyzed. All these patients underwent routine mapping and ablation by conventional approach [Non-Swartz sheath support (NS) approach] firstly. In the situation of the distal portion of the coronary venous system (CVS) not being accessed or a good target site not being obtained, the Swartz sheath support (SS) approach was attempted alternatively. If this still failed, the hydrophilic coated guidewire and left coronary angiographic catheter-guided deep engagement of Swartz sheath in GCV to support ablation catheter was performed. Results: A total of 103 VAs (103/123, 83.74%) were successfully eliminated in DGCV-AIV. By NS approach, the tip of the catheter reached DGCV in 39.84% VAs (49/123), reached target sites in 35.87% VAs (44/123), and achieved successful ablation in 30.89% VAs (38/123), which was significantly lower than by SS approach (88.61% (70/79), 84.81 % (67/79), and 75.95% (60/79), P < 0.05). There were no significant differences in complication occurrence between the NS approach and the SS approach (4/123, 3.25% vs. 7/79, 8.86%, p > 0.05). The angle between DGCV and AIV <83° indicated an inaccessible AIV by catheter tip with a predictive value of 94.5%. Width/height of coronary venous system>0.69 more favored a SS approach with a predictive value of 87%. Conclusion: For radiofrequency catheter ablation (RFCA) of VAs arising from DGCV-AIV, the SS approach facilitates the catheter tip to achieve target sites and contributes to a successful ablation.

Reconstructing gene regulatory networks in single-cell transcriptomic data analysis.

Gene regulatory networks play pivotal roles in our understanding of biological processes/mechanisms at the molecular level. Many studies have developed sample-specific or cell-type-specific gene regulatory networks from single-cell transcriptomic data based on a large amount of cell samples. Here, we review the state-of-the-art computational algorithms and describe various applications of gene regulatory networks in biological studies.